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"j" ' We use the LSZ reduction theorem and interpolating fields, alongwith 

the heavy quark effective theory, to investigate the structure of the Bethe- 

i-^ I Salpeter amplitude for heavy hadrons. We show how a simple form of this 

amplitude, used extensively in heavy hadron decay calculations, follows 



^' naturally upto 0(1/M) from these field theoretic considerations. 
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1 Introduction 

In recent years the heavy quark effective theory (HQET) |l[|, and the conse- 
quent heavy quark symmetry have emerged as useful tools in studying the decay 
properties of heavy hadrons. In our previous works - W6 have used the 
Bethe-Salpeter (B-S) formalism to derive the consequences of the heavy quark 
symmetry for weak transitions of hadrons of arbitrary spin. 

In the heavy quark mass limit we argued [§-0 that the Bethe-Salpeter (B-S) 
amplitude ( wave function ) for an arbitrary heavy meson, of momentum P and 
mass M, can be written in momentum space as 

where Xa^ is a projection operator which projects out a particular spin, parity 
state from an unknown orbital wave function A. Here pi and p2 are the momenta 
of the heavy and light quarks respectively and v is the four-velocity of the heavy 
meson defined as f = P/M. The momentum k is defined a.s k = pi — p2- 



The important observation is that x has the Bargmann-Wigner [p| - [12] form 



x{v^) = \{l + i)V{k). (2) 

with r(A;), a Dirac matrix, depending on the spin, parity and orbital angular mo- 
mentum 0. From this simple form of the B-S amplitude follow all the dramatic 
results about the reduction of form factors in heavy meson decays. 

We argued that this form of the B-S amplitude, eqs. (P and (P, arises 
because, in the heavy quark limit, in the leading order of the heavy quark effective 
theory (HQET) [|I| , [Q the heavy quark spin is decoupled from the light degrees of 
freedom and as a consequence the B-S amplitude satisfies the Bargmann-Wigner 
(Dirac) equation on the heavy Dirac index a. The baryons can also be represented 
in a similar form. 

In this short note we demonstrate that this form of the B-S wave function 
follows naturally, upto order 1/M, from the notions of interpolating fields and the 
LSZ reduction theorem, when combined with the zeroth order HQET. We will 
demonstrate this explicitly for the 0~ and 1~ heavy mesons and then show how 
it is easily generalised to arbitrary heavy meson resonances and heavy baryons. 
As an intermediate step we show that, in the heavy quark mass limit, only half 
the components of the quark fields appearing in the interpolating field contribute 
to the matrix element. For example, for a heavy {Qq) meson, in the rest frame 
of the meson only the quark part of the heavy quark field and only the antiquark 
part of the light quark field contribute. 

In Section ^ we discuss the consequences of the heavy quark mass limit for 
interpolating fields. In Section ^ we combine the results of Section |^ with the 
HQET to obtain the desired form of the B-S amplitude. 



2 Interpolating Fields and the Heavy Quark Mass 
Limit 

As we mentioned in the Introduction, a variety of interpolating fields can be used 
to represent a particular bound state within a reduction formula |l^. In this 
section we will first use the equations of motion to establish relations between 
certain classes of interpolating fields and the normalisation constants appearing 
in these interpolating fields. We will then look at these relations in the heavy 
quark limit. Then we will study the consequences of these relations when we use 
these interpolating fields in reduction theorems. As a first example consider the 
pseudoscalar heavy meson for which we can take as interpolating fields 

01 (x) = Zl^'^ —T4!g{x)^5^Q{x) (3) 

or 

Mx) = -ZT^-^Ti^,{x)^'l,M^) , (4) 

where the normalisation constants 

Ari = (0|T^,(0)75^Q(0)|P) (5) 

and 

Ar2 = -(0|T^,(0)^75^Q(0)|P) (6) 

are chosen to ensure the normalisation condition 

m,{x)\P) = Zl'\-^^-\ (7) 

for 2 = 1,2. 

One can now use the equations of motion 

{iP - m)ip = (8) 

for the quark fields to show that 

TUq +mqNi 



M N2 



and, further, from the normalisation condition eq. (|^ we get 

Ni M 



N2 ruQ + m, 



Q 



(9) 



(10) 



where uiq {niq) is the mass of the heavy (light) quark, leading to 0i = 02- Note 
that in the heavy quark limit, where mg + ruq -^ M, N1/N2 -^ 1. 



The pseudoscalar case is particularly simple and perhaps not very enlighten- 
ing. Much more interesting is the vector meson case which demonstrates some 
generic properties of interpolating fields which are not apparant in the pseu- 
doscalar case. We will, in fact, show that the corresponding two interpolating 
fields in the vector, 1~, case are not exactly equal but differ by a term which is 
of 0(-g) and can be ignored in the infinite mass limit. It will be then easy to 
show that such a situation always occurs for all meson resonances. 

An interpolating field (t)^{x) for a vector particle, with momentum P, should 
satisfy 

{Q\ct^,{x)\P) = Zl'\,e-^P-\ (11) 

where the polarisation vector e^ satisfies the transversality condition f ■ e = and 
e*^e^ = — 1. Thus, as in the pseudoscalar case, we can consider the following two 
interpolating fields 

4>l{x) = zl'^yTi,,{x)^^ijQ{x) (12) 

and 

0J(^) = -ZI''^-^T%{x)^'^^^^q{x) (13) 

where "la = l^i — ']^v^^. The normalisation constants, Vi, with i = 1,2, are given 
by 

V^ = -{0\MO)fMO)\P) (14) 

and 

Now one can use the equations of motion eq. (||) for the quark fields to obtain 

id^^g{x)Yl'l;^Qix) = -(mg + mg)[V'g(x)7^VQ(a;) - 2i^g{x) ^ ^qix)] ■ 

(16) 

— * — * — * 

where the transverse covariant derivative, D-~ = D^ — v^v ■ D, acts to the right. 
One sees that in this case the two interpolating fields are not exactly equal 
but differ by a term which can be neglected, as it is small in the heavy quark 
limit. Note that the term neglected, D-^/M, (as mq + niq -^ M) is exactly 
the kind of term which is dropped in deriving the HQET at leading order 0. 
Physically this means that we consider all momenta transverse to the line of 
flight of the meson to be small. One can also look at this result in another way. 

iil)q{x)j^ipQ{x) is in fact another possible candidate for the interpolating fleld 
for the vector meson but its overlap with the physical state becomes negligible 

in the heavy mass limit. Note also that iipg{x)^ipQ{x) corresponds to a p-wave 
contribution to the vector, 1~, state coming from the anti-quark part of the heavy 
quark. In the zeroth order HQET recall that the heavy quark and antiquark are 
described by separate flelds or to put it another way, the antiquark part of the 



heavy quark field is suppressed and thus the overlap of this interpolating field 
with the physical heavy meson state becomes very small in this limit. 

Hence, only in the heavy quark limit, with mg + niq — »• M, the two interpo- 
lating fields, (pj^ and 0^ are equal along with the important relation 

V^i = V^2, (17) 

between the normalisation constants. 

The above considerations turn out to be a general feature of the heavy quark 
limit. Given a particular interpolating field NiTtpqTtpQ then one can show that 
this is equal to the interpolating field N2^j^Tipg'y'^TiljQ, upto terms oiO{D-^/M), 
alongwith the relation between the normalisation constant Ni = N2, in the heavy 
quark mass limit. Here T is some Dirac matrix, possibly with derivatives. We 
shall come back to this point at the end of the next section. 

We now discuss the consequences of the equality of interpolating fields for 
matrix elements, in particular the Bethe-Salpeter amplitude, by using them in 
reduction theorems. First of all we would like to address the question : what 
does one mean that one is free to choose equivalent interpolating fields? We will 
address this question in terms of the Bethe-Salpeter amplitude but the conclusions 
apply to any matrix element involving a heavy meson or baryon. 

As an example we consider the heavy vector meson. We define the Bethe- 
Salpeter amplitude for the heavy vector meson as 

M^^{XuX2) = WQaiXi)GiXi,X2Wqix2)\P) , (18) 

In the above equations ipQ (V'g) is the heavy (light) quark field and \P) represents 
the heavy vector meson state with a certain momentum P and mass M. Here 
time-ordering is implicit. G{xi,X2) is a colour matrix chosen to make the B-S 
amplitude gauge invariant. However the following analysis also applies to non 
gauge invariant amplitudes. 

One can now use the reduction theorem to write the B-S amplitude as 

M,^(xi,X2) = lim fP^-M^) f d%e-'''-^0\ijQa{xi)G{xi,X2WJx2)(pHy)\0) . 

(19) 
The 0^ in ([l9|) is e'^^jj where 0^ is given in eq. (0) or eq. ([I3[) . Thus one can 
write 

Mj{x,^X2)=Yu.Zr^^^^- 

dV-'''-^(0|T^Qa(xi)G(xi,X2)^,^(x2)x/(a,)V'^(?/)^,,(|/)|0), 

(20) 

id 

where x{9y) is either /or -^^- Correspondingly, V in eq. (1^) is either Vi or V2. 



The statement that one can use any one of the alternative interpolating fields 
means that they should give the same answers in the reduction theorems. To see 

what this means consider the second interpolating field, i.e. x{dy) = -jf/in the 
matrix element (^). Integrating by parts we can shift the y derivative in x onto 
exp{—iP ■ y) to obtain 

M„^(xi,X2)= lim Zl^- ^■ 



d'ye-'^-y{0\Ti;Q^ix,)Gix,,X2W,{x2)ii//),'rQ{y)^Aym. 

(21) 

With the other choice of the interpolating field we have instead 



M« ^(xi, X2) = Jim^^ Z,' 

Ae"^^-^(0|r^Q„(xi)G(xi,X2)V^f(x2)(^/^^(l/)^,5(t/)|0). 

(22) 



Thus comparing we see, that if the two expressions, ( pTD and (^21) , are to be 
equal, then js/ = ^, as an operator inside the matrix element, in the mass shell 
limit. In general this is not much of a restriction, as Vi and V2 are functions 
of V. However, we have shown that, in the heavy quark limit, the ratio of the 
normalisation constants is unity leading to ^ = 1. In fact we see from ([21|) and 



(^21) that this condition means that only certain components of the quark fields 
entering in the interpolating fields contribute to the matrix element, i.e those 
satisfying 

ijQ{y)i = i^Q{y) 
My) = -^{y)- (23) 

within the reduction formula, in the mass shell limit. In other words, in the rest 
frame of the on-shell meson, only the quark part of the heavy quark field and 
surprisingly also only the antiquark part of the light quark field contributes to 
the matrix element. The same relation holds for the pseudoscalar case. Also 
the same relation will hold in any matrix element involving a heavy meson state 
when we use the interpolating fields to reduce the state. 

It is obvious that this physical requirement will be enforced by taking the 
interpolating field for the heavy vector meson, in the heavy quark limit, to be 

4>,{^) = ^(0i + 0j) 



In fact even if we start with an arbitrary linear combination of the two fields, 
only the sum survives in the heavy quark limit. 

With this choice the B-S amplitude picks up the correct projection operator 



1+^ 



to enforce the above condition and can be written in the form 



with 



and 



Mj{xuX2)=X,\v)Al^,{xuX2) 



X,\v) = -[{1 + ^)^,' 



(25) 



(26) 



a5\^li^2) — „,iii\^, ^3 



lim Z- 



,/2{P'-M') 



Vi 



Ae-'^-^(0|T^Q«(xi)G(xi,X2)^f(x2)^^(l/)^,5(y)|0) 



(27) 



The B-S amplitude for the heavy pseudoscalar particle is also of the above form 
with / replaced by 75 and Vi by A^i R. It is not very useful as it stands because 
of the unknown matrix function A^^. However in the next section we see that 
when combined with the HQET it leads to the previously proposed [^ - form 
of the heavy meson B-S amplitude. 



3 Bethe-Salpeter Amplitudes and HQET 

As is well known by now 0], one can go to the zeroth order HQET by simply 
replacing the QCD heavy fields in eq. ( P7| ) by the zeroth order heavy quark 
effective fields. Thus we shall now simply consider the heavy quark fields, ipq, 
appearing above to be the corresponding fields appearing in the effective theory. 
One can now decouple these heavy fields through the transformations H] 



W 



and 



where 



^q(x) = Q{x)W 



Q{x) 



W 



Pexp 



W 



dsA ■ V 



(28) 

(29) 
(30) 



is a path-ordered exponential, wherein the path is a straight line from — 00 to x 
along the w-direction 

^Although this form of the B-S amphtude has been forced on us in the heavy quark hmit, we 
could as well use it for the light meson as after all we are free to choose any linear combination 
of the two fields with the appropriate normalisation. 



Because the transformed fields, Q and Q are now decoupled from the gluons, 



the right hand side of eq. (^) factorises and we can now use the heavy quark 
propagator 



(o|ga(xi)Q^(i/)|o)= jd^p 

to write the B-S amplitude as 



Q-ip-{x-i-y) 
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(31) 



M„f^(xi,X2)= lim Z^^^ 



(P2 _ M' 



dVV-*''-'(0|^f(x2)Vr M G{X,,X2)W 



^~ip-{xi-y) 

]^v -p-niQ 



X,\v)i:gs{ym . 



(32) 



The projection operator x(f ) now gets rid of the ]^ in the heavy quark prop- 
agator to give 

Maf'{xi,X2) = Xa\v)Asf'ixi,X2) , (33) 



with 



A/(xi,X2)= \imZy 



1/2 (P2 - M2 



P2^M2 



/ d^yd 



^pe-iP-y. 



-ip-{xi-y) 

V ■ p-niQ 



Vi 



m^,{x2)W l^ G{x,,X2)W I " %,{ym 



(34) 



Eq. (^3[) is essentially the result we have been looking for. We can now 
transform this result to momentum space by defining the Fourier transform 



with 



Mj{p^,p2) = x{v)jA,\p^,P2) 



As^ipi,P2) = J d^xid*X2e'P'-'''e'P''''As^{xuX2) 



(35) 



(36) 



The B-S amplitude in this form now satisfies the Bargmann-Wigner equation on 
the 'heavy' label a, leading to all the nice results of the heavy quark symmetry. 
Further, from translation invariance one can write 

{0\i>^g{x2)W [l^] G{x^,X2)W [lY'^AvM = J d''qd^ke-"^<'''-^'^e-"'<^^-y^As''{q,k) 

(37) 
leading to 



As (pi,P2) = liiii Z-:> 



1/2 (-P^ ~ ^^)x4/D ^ / ^4, 

-b [P - pi - P2) d k 



As^{p2,k) 
V ■ {P - k) -rriQ 



{3i 



One way to generate the physical pole at P^ = M^ is for the Green's function 
A{p2,k) to be peaked at f ■ A; = A, where A = (M — rriQ). In other words 
the "longitudinal" mass of the light degrees of freedom is A. This is of course 
trivially true in the weak binding or equal velocity approximation, where there 
is no transverse momentum and m^ = A. 

It is now easy to generalise to arbitrary orbital resonances. One can take as 
interpolating fields, for orbital angular momentum L, either 

0i(^) = zl^'^T^p,ix)T^PQix) (39) 

or 

0i(x) = ~zl^'^^'-^T4^,{x)YrM^) (40) 

where the T are listed in the table for the four possible spin-parity states. We 
have omitted Lorentz indices and the Lj are normalisation constants. Again 
using the equations of motion it is easy to show that the difference between 

these two possible interpolating fields is always either of the form ipg- ■ ■ ^ .^ i^g 
or i^q' ' ' m +m "^Q ■ ^^ usual these are small in the heavy quark limit. (Here 

(9^= d^— d/i-) Hence in the heavy quark mass limit we find 0i(a;) = 0|(x) and 
Li = L2. Thus the rest of the analysis follows in general. One can see that it is 
also, in principle, quite easy to generalise to heavy baryons. 
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Table Caption. 

Lowest gamma structure ^ appearing in interpolating fields for mesonic L-wave 
states. 



"state J^^ 



'L,_i (L-i)w"^H-)"^^ ^ d:^...d 



M2 ^M-L 

± 
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